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Frustration in quantum many body systems is quantified by the degree of incompatibility between the lo¬ 
cal and global orders associated, respectively, to the ground states of the local interaction terms and the global 
ground state of the total many-body Hamiltonian. This universal measure is bounded from below by the ground- 
state bipartite block entanglement. For many-body Hamiltonians that are sums of two-body interaction terms, a 
further inequality relates quantum frustration to the pairwise entanglement between the constituents of the local 
interaction terms. This additional bound is a consequence of the limits imposed by monogamy on entanglement 
shareability. We investigate the behavior of local pair frustration in quantum spin models with competing inter¬ 
actions on different length scales and show that valence bond solids associated to exact ground-state dimerization 
correspond to a transition from generic frustration, i.e. geometric, common to classical and quantum systems 
alike, to genuine quantum frustration, i.e. solely due to the non-commutativity of the different local interaction 
terms. We discuss how such frustration transitions separating genuinely quantum orders from classical-like ones 
are detected by observable quantities such as the static structure factor and the interferometric visibility. 

PACS numbers: 03.65.Ud, 03.67.Mn, 05.30.Rt, 75.10.Pq 


I. INTRODUCTION 

Many body systems are typically modeled by global Hamil¬ 
tonians that are sums of local interaction terms. Frustration 
characterizes those situations for which it is impossible to sat¬ 
isfy simultaneously the minimization of all local interaction 
energies in the presence of the global constraint that imposes 
minimization of the total many-body energy M- 

In interacting many-body systems frustration may be due 
either to complex geometrical configurations and/or compet¬ 
ing interactions. These types of frustration are generic in 
that they are common to both classical and quantum sys¬ 
tems alike. On the other hand quantum systems, due to the 
non-commutativity of different local Hamiltonians, may fea¬ 
ture further, genuine quantum, sources of frustration with no 
classical counterpart I8 h1QI1 . In particular, quantum counter¬ 
parts of classically unfrustrated systems can be frustrated [§- 
i SMI. Therefore, for many-body quantum systems, a 
universal measure of total frustration encompassing all possi¬ 
ble sources, generic and specifically quantum, should quantify 
the degree of incompatibility between global and local ground 
states associated, respectively, to the total many-body Hamil¬ 
tonian and to the local interaction terms. 

Such a measure has recently been introduced, in Refs. Cl 
SI , exactly in terms of the overlap between the ground states 
of local interactions and the local reduced states obtained by 
partial trace from the global ground state of the total many- 
body Hamiltonian. It provides a unified treatment of both 
generic and specifically quantum contributions to the total 
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frustration of local interaction terms. Quantum frustration is 
then singled out by the existence of a quantum lower bound to 
the total frustration that is realized in terms of the ground-state 
entanglement between the local interacting subsystem and the 
remainder of the total many-body system. 

In the present work we first show how frustration relates 
not only to the global block entanglement across a biparti¬ 
tion of the total many-body system but also to the local entan¬ 
glement among the constituents of the individual interaction 
terms. Next, we show how these relations rule quantitatively 
the onset of ground-state dimerization and the formation of 
valence bond solids. For many-body Hamiltonians sums of 
two-body terms, the local interacting pairs are naturally iden¬ 
tified as the elementary subsystems. We deduce a direct re¬ 
lation between quantum frustration of local pairs and the pair 
entanglement as measured by the concurrence fl6l[l7ll . Here¬ 
with we are able to establish a relation between quantum frus¬ 
tration and entanglement monogamy lfl~8lfl~9h which imposes a 
strong constraint, without classical counterpart, on the share- 
ability of quantum correlations. 

Specifically, we investigate the interplay between generic 
and quantum frustration for an ample class of spin-1/2 mod¬ 
els with nearest-neighbor (NN) and next-to-nearest-neighbor 
(NNN) interactions. These models feature both quantum frus¬ 
tration due to the non-commutativity of the different local in¬ 
teraction terms and generic frustration due to competing inter¬ 
actions on different length scales. We show that ground-state 
dimerization in these systems is realized exactly at the tran¬ 
sition from generic to quantum frustration, i.e. at values of 
the Hamiltonian parameters for which the total frustration be¬ 
comes entirely due to the ground-state block entanglement. 

Remarkably, such models are amenable to quantum simu¬ 
lation with systems of trapped ions I2(il - l24ll . ultracold atoms 
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in optical lattices l25l - 281. coupled cavity arrays ||29l ,l30 : l, su¬ 
perconducting qubits IBlll32 ;l, or nuclear magnetic resonance 
simulators 1133113411 . Therefore, as we will show in the follow¬ 
ing, changes in the ground-state patterns caused by changes 
in the frustration properties can be experimentally observed, 
e.g. in optical-lattice realizations, by measuring the quasi mo¬ 
mentum distribution as quantified by the static structure fac¬ 
tor or the interferometric visibility lf37h . 

The paper is organized as follows. In Section [TT] we intro¬ 
duce the universal measure of total frustration in terms of the 
global-to-local infidelity and describe its main features. In 
Section [Til] we consider generic spin-1/2 models with pair¬ 
wise interactions on different length scales, discuss their frus¬ 
tration properties according to the degeneracy of the local 
ground spaces, establish rigorous conditions for the presence 
of quantum frustration, and derive an exact inequality on pair¬ 
wise frustration based on entanglement monogamy. In Sec¬ 
tion [IV] we use these results to investigate the transition from 
generic frustration to genuine quantum frustration in models 
with competing nearest-neighbor (NN) and next-to-nearest- 
neighbor (NNN) interactions, and identify different phases as¬ 
sociated, respectively, to generic frustration and to genuine 
quantum frustration. In Section [V] we discuss the observable 
characterization of such frustration transitions in terms of the 
static structure factor and the direct observation of the univer¬ 
sal measure of total frustration in terms of the interferometric 
visibility. Summary and outlook are discussed in Section lVTl 


II. FRUSTRATION 


Here we review the main aspects of a universal measure 
of total frustration expressed in terms of the incompatibility 
between local and global orders, and briefly discuss how this 
universal measure (global-to-local incompatibility, or global- 
to-local infidelity) connects to quantum entanglement. The 
class of many-body Hamiltonians Hr, under investigation 
comprises all those that can be written as sums of different 
local terms: Hq = he, where G stands for the global (to¬ 
tal) system and £ C G are the local subsystems associated 
to the local interactions he. Let us consider one such local 
subsystem £. Denoting the remainder of the total system by 
R = G\£, frustration of £ occurs when the local projection of 
the global ground state of Hq, obtained by partial trace over 
R, i.e. the reduced state of subsystem £, fails to belong to the 
local ground state space of he. Consequently, such frustration 
is directly quantified in terms of the overlap between the pro¬ 
jector TL e onto the local ground state space of the local Hamil¬ 
tonian he and the reduced local density matrix pe = Tt r {pg} 
from the global ground state pa of the total Hamiltonian 7T G . 
The ensuing universal measure of total frustration fe is thus 
defined as follows 0: 

// = 1 - Tr {(nr ® 1 R ) Pg } = 1 - Tr {U e pe} . (1) 

The measure is uniquely defined as long as G features a non 
degenerate global ground state space, while all local subsys¬ 
tems £ may feature local ground spaces of arbitrary degener¬ 
acy. 


If the global ground state is degenerate, then the measure 
can in general depend on which global ground state is actu¬ 
ally considered. Uniqueness is guaranteed by considering the 
a priori equiprobable statistical average, i.e. the convex com¬ 
bination of all possible degenerate global ground states with 
equal weights or, in other words, the maximally mixed global 
ground state (MMGGS) iflftl . 

Ground-state degeneracy in the presence of ordered phases 
with nonvanishing local order parameters is tied to sponta¬ 
neous breaking of some symmetries of the total Hamiltonian, 
while, by definition, the MMGGS preserves all symmetries of 
the model. Choosing a particular global ground state among 
the degenerate ones or the MMGGS depends on whether one 
is interested in features of the frustration that are, respectively, 
either state- or model-dependent. In the following we will 
study frustration in different quantum spin systems featuring 
both non-degenerate and degenerate global ground states. 

Majorization theory and the Cauchy interlacing theorem 
provide an immediate link between frustration and entangle¬ 
ment. Indeed, it is straightforward to show that fe is bounded 

( r ) 

from below by a quantity that depends uniquely on the 
reduced local density matrix pe 11 4 . 1 1 511 : 

ft > 4 r) , (2) 

where 

4 r) =i-i>i. (3) 

i —1 

The index r denotes the rank of the local ground space pro¬ 
jector lie and A' are the the first r largest eigenvalues of the 
reduced local density matrix pe in descending order. 

In the case of a non-degenerate global ground state pc the 

(r) 

quantity coincides with the Hilbert-Schmidt distance be¬ 
tween pc and the set of pure states with Schmidt rank less or 
equal to r 0. 

This quantity is a proper entanglement monotone | [38fl . 
Moreover, for r = 1, i.e. in the case of a non-degenerate 
ground space of the local interaction term he, it coincides with 
the bipartite geometric entanglement introduced in Ref. [[39]. 
This is the ground-state entanglement between the local block 
£ and the remainder R of the total many body system G as 
measured by the Hilbert-Schmidt distance Dhs of the pure 
ground state pc = |d/) G (4 , | G from the set of the pure bi- 
separable states |H/) bs defined as | v k) {ls = \ip)e ® \4>) R , so that 

4 0 = {| “i* } Dhs (I*>g > l*>6.) • 

Therefore 1 realizes a quantum lower bound to total frus¬ 
tration in terms of the bipartite geometric entanglement. 

The fact that the above measure of frustration is bounded 
from below by quantum entanglement when the many-body 
system features a pure, non-degenerate global ground state 
is a direct consequence of the fact that for pure states all 
quantum correlations reduce to entanglement. On the other 
hand, mixed states allow for quantum correlations more gen¬ 
eral than entanglement, such as the quantum discord, that can 
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be present even in separable states. One needs to keep track 
of this important difference when considering the general case 
of degenerate global ground states. Indeed, in order to intro¬ 
duce a measure of frustration that should not depend on the 
choice of one particular state among the possible degenerate 
ones, we need to introduce a reference global ground state that 
properly averages over all possible degenerate ground states 
and preserves all symmetries of the many-body Hamiltonian. 

This can be achieved either by taking a coherent quantum 
superposition of all the degenerate ground states with equal 
amplitudes, or by introducing the maximally mixed global 
ground state (MMGGS), i.e. the convex combination of all 
the possible pure ground states with equal probability weights. 
In such a mixed state, correlations are not completely charac- 

if) 

terized by quantum entanglement alone, and the quantity e) 
turns into the sum of the ground-state block entanglement be¬ 
tween £ and R, as measured by the convex roof of the pure- 
state bipartite geometric entanglement, plus the classical cor¬ 
relations that are established between the local subsystem £ 
and an ancillary system when the latter performs a general¬ 
ized quantum measurement on £ fl5l[40l l. 

Relation ([2} allows for a more refined classification of frus¬ 
tration in the quantum domain than merely distinguishing be¬ 
tween frustrated and non-frustrated systems, as in the classical 
domain. Namely, we can now introduce the following classi¬ 
fication of quantum systems according to frustration: 

1. Frustration-free systems (FF): 

f t = e (r) = 0 V f, 

2. Inequality saturating systems (INES): 

h = 4 ) > o vf, 

3. Inequality non-saturating systems (non-INES): 

ft > e^ r) > 0 Vf. 

Accordingly, if generic frustration is absent, the system is 
always in an INES ground state and thus either frustration-free 
if the ground state is classically correlated or genuinely quan¬ 
tum frustrated if it is quantum correlated. When the inequal¬ 
ity is not saturated (non-INES) the system’s total frustration 
in general comprises both generic and genuine quantum con¬ 
tributions. Indeed, as we will show in the following, the inter¬ 
play between generic frustration and genuine quantum frustra¬ 
tion and the transition from non-INES to INES configurations 
correspond to transitions from phases with classical-like mag¬ 
netic order to phases with genuine quantum order. 


III. FRUSTRATION AND MONOGAMY OF 
ENTANGLEMENT 

Let us consider translationally-invariant (periodic boundary 
conditions) many-body system consisting of N spin-1 /2 el¬ 
ementary constituents, featuring generic two-body (pairwise) 
interaction terms with the only constraint that they preserve 


the parity symmetry along the three independent spin direc¬ 
tions, so that the total many-body Hamiltonian reads 

u -»:>>?-XX. (5) 

i,j,ot ij 

where of are the Pauli spin operators ( a = x. y. z) for the 
i-th spin and J“ • is the interaction coupling strength between 
spins i, j along the a direction. For such Hamiltonians the ele¬ 
mentary local interaction subsystems £ are naturally identified 
as the ones made up by pairs of directly interacting spins with 
local Hamiltonian hij whose energy eigenstates are the four 
maximally entangled Bell states (singlet and triplet states). 

For such Hamiltonians that include, among others, the 
Ising, XY, and Heisenberg models, it is possible to introduce 
a quantum generalization of the classical Toulouse criteria for 
frustration-free systems mm. The first step of this con¬ 
struction consists in defining a prototype model. A quantum 
spin Hamiltonian of the type Eq. 0 is a prototype model if 
there exists at least one local ground state common to all lo¬ 
cal interaction terms hij and all pairwise couplings are ferro¬ 
magnetic (< 0). Having defined the prototype model, a 
quantum version of the Toulouse criteria can be formulated as 
follows: 

i) - The MMGGS of a prototype models is INES for all pair¬ 
wise interaction terms hij. 

ii) - For all models obtained from a prototype model by local 
unitary operations and partial transpositions acting on any 
arbitrary subsystem, the MMGGS is still INES on all pairwise 
interaction terms hij. 

The quantum Toulouse criteria allow to make predictions 
on the frustration of a model without requiring the complete 
and exact diagonalization of the system. In fact, they allow 
to state powerful exact relations between the degeneracy of 
the local ground spaces, entanglement, and the presence of 
quantum frustration in systems described by Hamiltonians of 
the form Eq. (0. Indeed, a key feature of a prototype model is 
that all interacting pairs admit at least one local ground state 
in common. Then, the following holds: 

Theorem. If a many body system with pairwise local interac¬ 
tions satisfies the quantum Toulouse criteria and the degen¬ 
eracy of the ground space common to all local interacting 
pairs has dimension r > 2, then the system is necessarily 
frustration-free, i.e.: ft = e'p =0 V£. 

Proof. If a model satisfies the quantum Toulouse criteria and 
the common ground space has dimension r = 2, then all in¬ 
teracting pairs admit as common ground state two out of the 
four maximally entangled Bell states. Obviously any linear 
coherent superposition of these two Bell states is again a legit¬ 
imate ground state of each local interacting pair. All possible 
superpositions include separable ones of the form |«) <g> |j). 
Hence, the fully factorized state | 'll fact) = ®\ = \ |*) is a le¬ 
gitimate ground state of the total many-body Hamiltonian, in 
close analogy to the situation that occurs at the factorization 
points of many-body Hamiltonians in transverse fields lf4ll — 
14311 . By immediate inspection, the ground state | ^ fact) has 
vanishing frustration, single-site entanglement, and pair block 
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entanglement for all interacting pairs. Moreover, since the to¬ 
tal energy associated to |4 , / ac t) is the sum of the minimum 
energies of all local pair interaction terms, any possible other 
total ground state must feature the same property and hence its 
projection onto any pair of interacting spins must belong to the 
pair local ground space, implying that its frustration vanishes 
on all pairs. Exactly the same argument applies also to the 
cases of local ground space degeneracy r = 3 and r = 4. □ 

We are thus left with the only non-trivial case of a local 
ground space degeneracy r = 1, corresponding to local pairs 
of interacting spins admitting one of the four Bell states as a 
local ground state. In such a case, the quantum lower bound to 
frustration always coincides with the ground-state geometric 
entanglement between the local block (pair) £ = { i. j } and 
the rest R= G\£ of the system. Dropping from now on the 
index r = 1, it reads 


€£ — 1 A max ; (6) 

where \ max is the largest eigenvalue of the reduced local 
density matrix pi. From now on, without loss of general¬ 
ity, we will consider the case of anti-ferromagnetically in¬ 
teracting pairs of spins admitting as unique, non-degenerate 
ground state the antisymmetric singlet state: I?//- 5 ) = 

-^= (| t4j) — | 4-itj})- The corresponding projector onto the 
local ground space then reads 

= \tp ij ){tp ij \. ( 7 ) 

Exploiting the symmetry under parity along the three spin di¬ 
rections, the two-body reduced local density matrix takes the 
general form 

pe=ij = lij + Y Sij a i a j • ( 8 ) 

a 

Here gf denote the spin-spin correlation functions: gf = 
(SfSf), where Sf = (1/2) are the spin operators and 
(•) indicates the expectation value over the ground state of the 
system. Consequently, the frustration of each interacting local 
spin pair reads 

/« = ! + £»«• w 

ot 

and the quantum lower bound e l3 on the total frustration ])., 
of each local interacting pair reads 

3 

<rj = -y ' min { !f'ii + tii, + 9ij,9fj - 9ij ~ 9?j , 

~9ij + 9ij - 9ij , ~9ij ~ 9ij + 9ij } ■ (10) 

Since the eigenvectors of pt—ij are the four Bell states we can 
establish a direct connection between the total frustration f l3 , 
the quantum lower bound e,j, i.e. the ground-state entangle¬ 
ment between the local interacting pair ij and the remainder 
of the systems, and the entanglement of formation between 
the spins i and j in the pair, as measured by the concurrence 
Cij lfl6LrL7il . We have that the following inequality holds: 

Cij = max(0,1 - 2eij) > max(0,1 - 2 fa) , ( 11 ) 


where Cij is the concurrence of the spin pair {/, j} and we 
have exploited the explicit analytic form that it takes on the 
set of the four Bell states. Equality holds only when the frus¬ 
tration is INES, i.e. when . Therefore, saturation of 

the inequality implies a direct correspondence between quan¬ 
tum frustration, a global feature of the pair as measured by its 
block geometric entanglement with respect to the remainder 
of the entire many-body system, and the local quantum entan¬ 
glement within the pair. The more frustrated the subsystem 
{i,j} is with respect to the remainder of the many-body sys¬ 
tem, the weaker is the internal entanglement between spins i 
and j forming the pair. Moreover, since the quantum lower 
bound ejj is the entanglement between the local pair and the 
remainder of the many-body system if and only if the global 
ground state is pure, from Ineq. (ITU it follows for all pure 
global ground states that the stronger the entanglement be¬ 
tween spins in the local pair, the weaker is the entanglement 
of the pair with the remainder of the many-body system. 

These relations allow for a restatement of frustration in 
terms of entanglement shareability. Indeed, exploiting the 
Coffman-Kundu-Wootters-Osborne-Verstraete monogamy in¬ 
equality lfl8[ IT9h that sets precise constraints on the entan¬ 
glement that can be shared among many spins, we find the 
following “frustration monogamy”: 

Y ( max (°’ 1 “ 2 fa )) 2 < Y C h 

j&Li j£Li 

< Y C lo < n ■ d2) 

jeG\i 


In the above, L, denotes the subset of spins that interact with 
spin i, while G\i denotes the total spin system without spin 
i. The upper bound r, is the tangle Iil8l 1 19( 1. When the global 
ground state is non degenerate, it coincides with the squared 
concurrence measuring the single-site entanglement of spin i 
with all other spins in the system. For spin systems described 
by Hamiltonians of the form Eq. ([5]) the tangle takes the ex¬ 
pression Tj = 1 — mf, — rriy — ml, where m a = (of) is the 
local magnetization along the a direction. Therefore r, = 1 
in the MMGGS, since the latter preserves all symmetries of 
the total many body Hamiltonian and therefore does not al¬ 
low for a spontaneous magnetization along any spin direction. 
This fact and the frustration monogamy Ineq. (Il2l) imply ex¬ 
act bounds on the minimum possible amount of frustration in 
a many-body system. Considering for simplicity the uniform 
case fj = f V{«, j}, and denoting by Nl the number of 
spins interacting with spin i, we have: 




(13) 


We thus find that the minimum frustration threshold varies 
from fmin = 0 for Nl = 1 for a many-body system con¬ 
sisting in just a single pair of interacting spins, to / min = 1/2 
at the thermodynamic limit of fully connected many-body sys¬ 
tems, i.e. with Nl —> oo. 
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FIG. 1: (Color online) Behavior of Ineq. 10 for the local pair frustra¬ 
tion as a function of the admissible values of the correlation functions 
gfj = g'lj and gfj in the reduced local pair density matrix pi = ij of 
an XXZ model with competing NN and NNN interactions. Verti¬ 
cally striped green area: reduced states pe=ij with non-INES frus¬ 
tration (strict inequality, generic frustration). Horizontally striped 
orange area: reduced states with INES frustration (inequality satu¬ 
rated, genuine quantum frustration). Red area: disentangled reduced 
states with INES frustration (inequality saturated, vanishing pair con¬ 
currence, enhanced genuine quantum frustration). Remaining white 
area: non-physical states. 


FIG. 2: (Color online) a) Schematic structure of a one-dimensional 
quantum spin model with competing NN (Ji) and NNN (J 2 ) 
anti-ferromagnetic interactions (periodic boundary conditions), b) 
Ground-state structure at the integrable Majumdar-Ghosh point 
J 2 /J 1 = 1/2. At this point the model features two exactly dimer¬ 
ized ground states that are the tensor product of NN spin singlets. 
The two valence bond solids break the translational invariance of the 
Hamiltonian with respect to the lattice spacing a. 


IV. FRUSTRATION AND DIMERIZATION 

The frustration f, :) of the local pair interaction //,, y is in gen¬ 
eral a function of the parameters entering in the total many- 
body Hamiltonian. Moreover, the interplay between different 
sources of frustration may lead to significant changes in the 
structure of the global many-body ground state. In particular, 
we wish to investigate the changes in the structure of quan¬ 
tum ground states when frustration turns from non-INES to 
INES, i.e. when genuine quantum frustration overrules com¬ 
pletely other generic sources of frustration. We will focus on 
the paradigmatic class of frustrated one-dimensional models 
with competing nearest-neighbor (NN) and next-to-nearest- 
neighbor (NNN) anti-ferromagnetic interactions, the so-called 
Ji — J 2 models 0: 

H = cosgbJ2 cos 6(SfS? +1 + S?S? +1 ) + sin 8S?S? +1 

i 

+ Scos<5 (Sf Sf +2 + Sf Sf +2 ) + sin SS?S? +2 .(14) 

i 

These systems are characterized by the ratio between NNN 
and NN interaction strengths J 2 /J 1 = tan</>, with (f> £ 
[0, |]), and the ratio between anisotropies tan S, with S £ 
[0, tj]. The translationally invariant case (periodic boundary 
conditions) is studied by exact diagonalization with N = 24 
spins. In Fig. Q] we report the behavior of the pair frustration 
fij , according to the classification induced by Ineq. (0. as a 
function of the spin-spin correlation functions, where we have 
exploited Eqs. ®, ©, and CEB, and the symmetry g = gf^. 
Both non-INES (generic) and INES (genuine quantum) frus¬ 
trations are realized in a significant range in the parameter 
space. 

In Fig. 0 we report schematically the structure of the quan¬ 
tum ground states for models with competing NN (Ji) and 


NNN (J 2 ) interactions, in the class of Eq. CLS. Exact ground- 
state dimerization is realized at the Majumdar-Ghosh point 
( t>MG — 0.45 solution of tan <f> = 1/2 [45]. This is the only in¬ 
tegrable point of models in the class Eq. (IT4l) . At this point the 
ground state becomes doubly degenerate, breaking explicitly 
the translational invariance of the lattice, which is recovered 
at double lattice spacing 2 a. The two ground states are tensor 
products of NN spin singlets, realizing the first and simplest 
paradigmatic instance of valence bond states ll46l - l48tl . Except 
for the Majumdar-Ghosh point, the global ground state is al¬ 
ways unique. At 4>mg frustration is evaluated in each of the 
two dimerized ground states as well as in the MMGGS, i.e. 
their convex combination with equal weights. 

In Figs.[3]and[4]we report the behavior of the pair frustration 
fij, the quantum lower bound f ., :) , and the pair concurrence 
Cij as functions of cj >, respectively for the NN (i,i + 1) and the 
NNN (i,i + 2) local pair interaction terms. The plots are taken 
at anisotropies 5 = 0 (XX symmetry), ^ (XXZ symmetry), 
and S = (classical Ising symmetry). The fully isotropic 
Heisenberg case (XXX) realized by S = j is reported in 
Fig.[5]for which we provide also insets comparing the detailed 
behavior of the pair frustrations in the MMGGS and in the 
two degenerate dimerized pure ground states at the Majumdar- 
Ghosh point 4>mg- 

Remarkably, from Figs.[3]and[4] we see that for all quantum 
models (XX, XXZ, and Heisenberg) it is always fi,i +1 = 
e* j+i, i.e. frustration of the NN interacting local pair is al¬ 
ways INES (genuine quantum), notwithstanding the presence 
of a source of frustration, due to geometry, in terms of com¬ 
peting interactions on different length scales. 

The NN pair frustration increases for increasing relative 
strength of the NNN interactions, reaching asymptotically the 
maximum value 3/4. This coincides with the maximum pos¬ 
sible value of e^j+i, i.e. the maximum possible value of the 
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ground-state block entanglement between the NN pair and the 
rest of the system. Such a maximum is achieved for a reduced 
pair density matrix pi^+i with all four eigenvalues degenerate 
and equal to 1/4. 

In other words, maximizing the relative strength of the 
NNN interaction maximizes the ground-state block entangle¬ 
ment e^i+i between NN pairs and the rest of the system 
and the corresponding (genuine quantum) NN pair frustration 
fi i+i- Hence, by monogamy, it minimizes the pairwise entan¬ 
glement Ci t i +1 between the spins forming the NN interacting 
local pair. Indeed, for a sufficiently large value of the ratio 
J2/J1, the NN pair concurrence C,.,+i vanishes exactly and 
remains zero thereafter. 

Even more important. Figs. [3] and Q] show that frustration 
fi,i+2 of the NNN interacting local pair undergoes a full tran¬ 
sition from non-INES (generic frustration) to INES (genuine 
quantum frustration) exactly at the Majumdar-Ghosh point 
tan (j> = 1/2. Therefore, transition to exact ground-state 
dimerization corresponds to exact transition from non-INES 
to perfect INES frustration of the NNN local interaction terms. 
For increasing relative strength of the NNN interactions the 
NNN pair frustration decreases, reaching asymptotically its 
minimum (whose actual value depends on 8). Such transitions 
in frustration exhibit universal features in the XX, XXZ, and 
Heisenberg cases, irrespective of the corresponding different 
symmetry classes. 

As shown in Figs. 0 and HI the frustration pattern differs 
radically only in the Ising case, as should be expected, since 
all local interactions commute in Ising systems, entanglement 
vanishes, and, as discussed in Section |III1 there are no gen¬ 
uine quantum sources of frustration. In this latter case we 
have simply two well distinguished regimes separated by the 
Majumdar-Ghosh point 4>mg- Indeed, at this point the NN 
frustration jumps from zero to half the maximum value 1/2 
and remains constant thereafter, while the NNN frustration 
jumps from the maximum value 1 to a vanishing value and 
remains zero thereafter. 

On the other hand, in the nontrivial quantum cases with 
non-commuting local interactions (XX, XXZ, Heisenberg) 
we observe a rich structure with three different patterns of 
frustration as the ratio of NNN to NN interactions varies. At 
small <fi well below the Majumdar-Ghosh point <Pmg (domi¬ 
nating NN interactions) NN spin pairs are weakly frustrated 
in opposition to NNN pairs that tend to be strongly frustrated. 
Indeed, for </> —0 the NNN pair frustration achieves values 
very close to unity, i.e. the maximum possible value, which is 
achieved exactly in the classical Ising case. For </’s far above 
the Majumdar-Ghosh point 4>mg (dominating NNN interac¬ 
tions) the situation reverses: NNN pairs become weakly frus¬ 
trated while NN pairs become strongly frustrated. 

These two limiting regimes are characterized by two simple 
antiferromagnetic orders that essentially coincide with those 
of the classical Ising case. Fig. [6] summarizes this behav¬ 
ior pictorially in terms of the lattice structure. In the limit 
<j> 0 the lattice reduces to a single chain with lattice spac¬ 

ing a and NN anti-ferromagnetic interactions, allowing for a 
simple anti-ferromagnetic order. In the limit 0 -0 tt/ 2 the lat¬ 
tice effectively splits into two chains, each with lattice spacing 



FIG. 3: (Color online) Frustration fi,i +1 of the NN interacting local 
pair, quantum lower bound eqi+i, and NN pair concurrence C/i+i, 
as functions of the parameter <f> ruling the ratio tan (j> between NN 
and NNN interactions. Upper panel: Ji — J 2 model with XX sym¬ 
metry, (5 = 0. Solid green line: /;,i+i = eqi+i, as the two quantities 
in this case coincide V0. Dashed blue line: C/i+i. Central panel: 
Ji — J 2 model with XXZ symmetry, 8 = n/6. Also in this case 
fi,i +1 = U,i+ 1 V0. Both are denoted by the solid green line. Dashed 
blue line: Ci,i+ 1 . Lower panel: Ji — J 2 model with Ising symmetry, 
8 = 7t/ 2. Dot-dashed red line: /qi+i between the Majumdar-Ghosh 
point (j>MG and 4> = 7t/2. Solid green line: fi.i+i = eqi+i = 0 
between the Majumdar-Ghosh point 4>mg and <f> = 0. Dashed blue 
line: Cqi+i. NN frustration is always INES. fi,i +1 = £i,;+i, for 
the quantum XX and XXZ models. In the regime of dominating 
NNN interactions, fi,i+ 1 —> 3/4, which is the maximum possible 
block entanglement eqi+i between the NN interacting local pair and 
the rest of the system. The leftmost vertical line corresponds to the 
Majumdar-Ghosh point 4>mg- The rightmost vertical lines corre¬ 
spond to the boundary value 4>b of the intermediate quantum regime: 
</>£, ~ 0.99 and (pt, ~ 1.08, respectively for the A'A' and the XXZ 
symmetries. Models with Ising symmetry feature only a frustration 
due to geometry. It undergoes a sharp transition from fi,i+i = 0 to 
fi,i +1 = 1/2 exactly at <Pmg and remains constant thereafter. The 
quantities e»,i +1 and C;,i +1 are always vanishing in the Ising case. 


2 a and NN anti-ferromagnetic interactions. Two simple anti¬ 
ferromagnetic orders are then established on each chain. 

However, in the quantum models a further intermediate 
regime is established in the central interval starting at the 
Majumdar-Ghosh point and extending to the right, i.e. with 
increasing values of </>. In this interval, we observe a smooth 
crossover in terms of a succession of small but sizable step- 
jumps in the behavior of the NN and NNN pair frustrations. 
This intermediate, central regime corresponds to the onset of 
a quantum order that is well distinguished from the two clas¬ 
sical (Ising-like) anti-ferromagnetic regimes discussed above. 
The lower and upper values of this interval are, respectively, 
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FIG. 4: (Color online) Same as in Fig.[3]but this time for the frustra¬ 
tion fi,i+ 2 , quantum lower bound ti,-i.+ 2 , and concurrence Ci,i +2 of 
the NNN interacting local pairs. Transition to genuine quantum frus¬ 
tration (non-INES to INES), fi,i+ 2 = e%,i+ 2 , occurs exactly at the 
Majumdar-Ghosh point. Frustration remains INES thereafter. In the 
regime of dominating NNN interactions, +2 —>■ 1/4, which is the 
minimum possible block entanglement ti,i +2 between the NNN in¬ 
teracting local pair and the rest of the system. Models with Ising 
symmetry feature only frustration due to geometry. It undergoes 
a sharp transition from fi,%+ 2 = 1 to fi,i+ 2 = 0 exactly at the 
Majumdar-Ghosh point and remains vanishing thereafter. 


the Majumdar-Ghosh point and the boundary point 4>h cor¬ 
responding the last step-jump in the increasing behavior of 
and decreasing behavior of fi.i+ 2 - 

The fact that a sharp change in the behavior in the frustra¬ 
tion properties appears exactly at the Majumdar-Ghosh point 
may seem at first surprising, since the transition to the dimer¬ 
ized phase in the J\ — J 2 model, regardless of the anisotropy, 
occurs at a critical value </> = <f> c much lower than (1>mg- For 
example, in the case of the fully isotropic Heisenberg model, 
S = 7t/ 4, the critical value <p c ~ tan _1 (0.2411) ~ 0.237, 
as determined by extensive numerical investigations based on 
the density matrix renormalization group (DMRG) and other 
numerical algorithms [54-4;]. However, the frustration does not 
feature any sharp change when crossing </> c . 

The smooth behavior of the frustration at <j> c and the appear¬ 
ance of step-jumps in the frustration when the system is in the 
intermediate quantum regime deserve further comments. Let 
us first consider the physical origin of the step-jumps. One 
might conjecture that they are due to finite-size effects and that 
they should disappear when considering much larger chains. 
In fact, this is not the case, as can be already understood by 
considering chains of relatively small size. Indeed, by looking 
at the structure of the ground state in correspondence to the 
occurrence of the step-jumps, we have verified by analyzing 


+ 

u 

t. 

+ 


+ 


+ 
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4>mg kH 

4> 


n/2 


FIG. 5: (Color online) Same as in Figs.[3]and Fig.[4]but this time for 
Ji — J 2 models with Heisenberg {XXX) symmetry, $ = 7r/4. In 
this case the boundary value of the intermediate quantum region is 
(j>b — 1.18. Insets: zoom of the behavior of the NN and NNN pair 
frustrations at and around the Majumdar-Ghosh point. Away from 
4>mg the pure global ground state is unique and thus always coin¬ 
cides with the MMGGS. Exactly at 4>mg, besides the averaged ones 
in the MMGGS, we also report the non-averaged NN and NNN pair 
frustrations in the two exactly dimerized pure global ground states. 
The NN pair frustration fi,i +1 vanishes in the ground state for which 
{*, i + 1} belong to the same dimer and are thus maximally entan¬ 
gled between themselves (purple full circle). Vice versa, it coincides 
with the block entanglement ti,i +1 at its maximum value 3/4 in the 
ground state for which {i, i + 1} belong to different dimers (green 
full square). The NNN pair frustration fi,i +2 = 3/4 (green full cir¬ 
cle) in both ground states, since {i, i + 2} always belong to different 
dimers. 


the numerical data obtained from exact diagonalization that 
at these points exact crossovers are realized between ground 
states of different parity. These crossovers are then respon¬ 
sible for the sudden increasing of the NN frustration and the 
corresponding sudden decreasing of the NNN frustration that 
are observed in the interval \4>mGi <t>b\- 

This structure of ground-state crossovers appears similar to 
the one occurring in XY chains in transverse field h, for val¬ 
ues of h below the factorizing field hf ifidi -Poll. Indeed, one 
can draw a very close parallel between the two-site factoriza¬ 
tion of the ground state of the J\ — J 2 model at the Majumdar- 
Ghosh exact dimerization point and the single-site fac¬ 

torization of the ground state of the XY model at the exact 
factorizing field hf. Firstly, in both cases, the sharp change in 
the behavior of various physical quantities, including frustra¬ 
tion and ground-state entanglement, that occurs exactly at the 
respective factorization points (single-site factorization in the 
XY model at hf, two-site factorization in the J\ — J 2 model 
at 4>mg)i is always present, independently of the size of the 
system [4311491 DO! . Moreover, in the finite-size XY chain, 
the crossovers occur between ground states of different par- 
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ity and appear for 0 < ft < hf. Such crossovers are always 
present, but their number and position depend on the size of 
the system 0. The crossovers terminate exactly at ft = 0 
and ft = hf. At each crossover, several physical quantities, 
including frustration and ground-state entanglement, feature 
step-jumps completely analogous to the ones that are featured 
in the central quantum regime of the finite-size — J -2 model, 
i.e. in the region 4>mg < 4 1 < <?V The analogy is complete by 
observing that the crossovers terminate exactly at <j> = (!>mg 
and (j> = (f>b, and that both exact ground state factorization 
and exact ground state dimerization are not associated to any 
quantum phase transition. Indeed, as already mentioned, both 
the Majumdar-Ghosh point in the J\ — J 2 model and the fac¬ 
torization point in the XY model occur independently of the 
size of the system. Moreover, both points are always associ¬ 
ated to a double degeneracy of the global ground state. 

In the XY model, where all quantities can be evaluated 
analytically and controlled precisely both at finite-size and in 
the thermodynamic limit, it is possible to verify exactly that 
the number of the step-jumps increases while the height of 
the jumps decreases as the size of the chain increases; in the 
thermodynamic limit the height of the jumps vanishes exactly 
and a fully analytic behavior is restored. Unfortunately, the 
non-solvability of the J\ — J 2 model prevents us, at least at 
the moment, to establish whether the central quantum region 
extends indefinitely, i.e. the boundary point increases as the 
size of the chain diverges, and whether analyticity is restored 
in the thermodynamic limit. 

Concerning the smooth behavior of the frustration observed 
at 4> c in small chains, we remark that the dramatic change to 
exact ground state dimerization occurs at and only at 4>mg- 
Furthermore, our investigation relies on exact diagonalization, 
which is feasible only for chains of rather limited size. To ad¬ 
dress larger sizes, one would need to resort to DMRG and 
related algorithms such as matrix product states (MPS). How¬ 
ever, the existing codes are not suitable, as they stand, for 
the evaluation of the frustration measure. The first stumbling 
block is that the DMRG algorithm needs profound modifica¬ 
tions and improvements in order to compute reduced density 
matrices and projectors onto the local ground spaces, their 
eigenvectors, and their spectra, in the presence of periodic 
boundary conditions. The second stumbling block is that, 
while DMRG allows to evaluate without much effort the cor¬ 
relation functions between neighboring sites, it works much 
less efficiently when one needs to compute correlation func¬ 
tions between non-neighboring spins in the chain. Unfortu¬ 
nately, exactly such correlation functions are those needed 
in order to evaluate, for instance, the static structure factor. 
These difficulties have motivated us, whenever direct exact 
analytical evaluation proved to be impossible, to resort to ex¬ 
act diagonalization algorithms based on augmented Lanczos 
methods. 

In case of sufficiently large chains, the frustration mea¬ 
sure is indeed expected to detect the transition to approxi¬ 
mate dimerization by showing an inflection at <j> c , in the same 
way as it does in exactly solvable models that we are cur¬ 
rently verifying, such as the 1-D XY model, or the 1-D Ki- 
taev and cluster-Ising models of symmetry-protected topolog- 



FIG. 6: (Color online) Lattices effectively realized in the different 
limits of cf>. For cj> 0 NNN interactions are strongly suppressed and 
the lattice reduces to a NN interacting chain with lattice spacing a. 
For <j> —¥ 7t/ 2 the lattice splits into two independent NN interacting 
chains, each with lattice spacing 2a. 


ical order 0. However, in the case of the J\ — J 2 model 
this behavior can become evident only once more powerful 
numerical techniques are developed in order to evaluate the 
frustration in chains of sizes much larger than the ones that 
are allowed by resorting to exact diagonalization. 


V. OBSERVABILITY: STRUCTURE FACTOR AND 
INTERFEROMETRIC VISIBILITY 


The nature of the intermediate, central quantum regime in 
comparison to the two standard anti-ferromagnetic orders is 
best captured by studying the behavior of the static structure 
factor, as defined below and reported in Fig. [7] Indeed, as 
changes in the patterns of quantum frustration imply a tran¬ 
sition between different lattice geometries, we can connect 
them to the quantum momentum distribution observed in ex¬ 
periments, which is quantified by the static structure factor 

Sf(k ): 


O JV 

s f( k ) = ^ cos M*-il)<^-^> > (!5) 

i,j =1 

where a is the lattice spacing, k is the wave vector, and S* is 
the spin operator vector of the ith whose components are the 
three Pauli spin-1/2 matrices. In most experimental situations 
direct access to the correlation functions for each pair can be 
extremely challenging. In such cases an “averaged” and “col¬ 
lective” information such as the static structure factor, that can 
be obtained from the analysis of the visibility which quantifies 
the contrast of the time-of-flight images 135113611 . may provide 
useful coarse-grained information on the changes in the quan¬ 
tum orders due to changes in the frustration patterns. 

In Fig. [7] we report the behavior of the static structure fac¬ 
tor for three different values of qb, respectively well below, 
around, and well above the Majumdar-Ghosh point 4>mg — 
0.45, comparing the quantum XX and Heisenberg models 
with the classical Ising case. For </> ~ 0.14, well below 
4>mg quantum and Ising models all features a single peak at 
k = it I a, corresponding to the same anti-ferromagnetic or¬ 
der. At <j> ~ 1.4, i.e. well above (f>MG a ll models feature two 
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FIG. 7: (Color online) Behavior of the quasi-momentum distribution 
Sf(k) as a function of the quasi-momentum k for different values 
of <j>. Black dot-dashed curve: classical Ji — J 2 model with Ising 
symmetry ( Sn/2 ). Red solid curve: quantum Ji — J 2 model with 
XX symmetry (8 = 0). Blue dashed curve: quantum Ji — J 2 
model with Heisenberg symmetry (5 = 7r/4). Upper panel: dom¬ 
inating NN interactions ( <j> ~ 0.14). The static structure factor 
features a single peak for all models, corresponding to a simple 
Ising-like anti-ferromagnetic order in a simple chain with NN anti¬ 
ferromagnetic interactions and lattice spacing a. Lower panel: dom¬ 
inating NNN interactions (</> ~ 1.4). The lattice splits in two sub¬ 
chains, each with lattice spacing 2a. Correspondingly, two symmet¬ 
ric peaks are featured by all models, corresponding to simple anti¬ 
ferromagnetic order on each sub-chain. Central panel: intermediate 
regime (<j) ~ 0.55). In this regime, the classical model still fea¬ 
tures two simple anti-ferromagnetic orders, while the quantum mod¬ 
els feature a nontrivial quantum order and no peaks (the many small 
secondary peaks are due to finite-size effects that vanish in the ther¬ 
modynamic limit). 


symmetric peaks at k = n /2a and k = 37r/2a, correspond¬ 
ing to two anti-ferromagnetic orders, one for each of the two 
effective chains of Fig. [6] These are the geometric regimes 
for which classical Ising and quantum XX and Heisenberg 
orders are indistinguishable. 

Finally, for <f> ~ 0.55, a value above but comparable to 
(f>MG, we observe that the quantum models do not feature any 
peak, corresponding to a quantum order totally distinguished 
from the two classical anti-ferromagnetic orders featured by 
the Ising model on the two effective chains. 

On the other hand, although the analysis of the static struc¬ 
ture factor provides useful information about the behavior of 
the frustration measure, in general it cannot yield complete 
information about the relation between the frustration and its 
quantum lower bound e^ . Indeed, while the existence of 
a direct relationship between the static structure factor and 
the frustration is due to a qualitatively similar dependence of 


both quantities on sums of correlation functions, a fact that 
can be immediately appreciated seen comparing Eq. © with 
Eq. (fl5] >. a direct connection between correlation functions 
and eij is hindered by the highly nonlinear nature of their mu¬ 
tual relationship, where the nonlinearity originates from the 
nontrivial minimization that occurs in Eq. ( ITOl i. 

As already mentioned, the genuine quantum regime estab¬ 
lished at the Majumdar-Ghosh point 0mg extends up to a crit¬ 
ical boundary value <j> b corresponding to the last step-jumps 
in the pattern of increasing NN pair frustration and decreasing 
NNN pair frustration, as reported in Figs. [3] [4] and [5] The 
extension of the quantum region is a non-universal feature, as 
it increases, as should be intuitively expected, with increasing 
symmetry of the different quantum models. Indeed, for <5 = 0 
(XX), (j) b ~ 0.99. For 5 = tt/6 (XXZ), cj> b ~ 1.08. Finally, 
for 6 = 7t/ 4 (Heisenberg), (j> b ~ 1.18. 

Concerning observability, besides measuring the static 
structure factor in time-of-flight experiments, a complemen¬ 
tary strategy would be the measurement of the visibility in 
interferometric experiments with atoms lf37tl - In quantum op¬ 
tics, the interferometric visibility V by feeding two successive 
photons from a source into a 50:50 beam splitter is 

V = Ti ( PaPb ) , (16) 

where p a and p b are the states of the two incoming photons. 

In fact, it is possible to envisage very general schemes for 
the estimation and the direct measurement of linear and non¬ 
linear functionals of quantum states, based on simple quan¬ 
tum networks ll52h . Such schemes have been recently in¬ 
vestigated and specialized to quantum simulators realized 
with atomic ensembles in optical lattices, either in terms of 
quantum switches |;53j] or multiparticle atomic interferome¬ 
ters [54], and are close to experimental realizability & 

VI. CONCLUSIONS AND OUTLOOK 

In summary, we have shown that a measure of frustration 
based on the global-to-local incompatibility overlap (global- 
to-local infidelity) quantifies efficiently both generic (geomet¬ 
ric, common to and equal for classical and quantum systems 
alike) and genuinely quantum sources of frustration. 

For a large class of Hamiltonians we have proven that the 
presence of genuine quantum frustration is possible only if the 
ground states of the local interaction terms are non degenerate, 
that is, the rank of the projection operator onto the ground 
space of the local interacting subsystems is r = 1. 

In such cases the total frustration can be related, via differ¬ 
ent upper and lower bounds, with two types of entanglement: 
the internal (local) entanglement between the constituents of 
the local interaction terms, and the block (global) entangle¬ 
ment between the local interaction terms and the remainder of 
the total many-body system. Moreover, we have established a 
“monogamy of frustration” relation quantifying the trade-off 
on frustration due to the interplay between local and global 
entanglement. 

We have then applied these general results to spin Hamilto¬ 
nians with competing anti-ferromagnetic interactions on dif- 
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ferent length scales, thus featuring simultaneously generic and 
genuine quantum contributions to the total frustration. We 
have then shown that the quantum phase transition to ex¬ 
act ground-state dimerization and valence bond solids corre¬ 
sponds to an exact transition to genuine quantum frustration. 

We have investigated the behavior of the total frustration 
as a function of the physical parameters of a generic class of 
many-body Hamiltonians. We have identified three different 
regimes, two generic ones in which frustration due to geome¬ 
try dominates and only one of the competing interaction terms 
is strongly frustrated while the remaining one is essentially 
frustration-free, giving rise to phases with a well defined, 
classical-like, anti-ferromagnetic order, and an intermediate 
one in which genuine quantum frustration is non-negligible, 
the total frustrations of the competing interaction terms are 
comparable, and a quantum valence bond solid phase is estab¬ 
lished with no definite magnetic order. 

We have further shown how these different many-body 
regimes can be clearly identified in an experiment by inves¬ 
tigating the behavior of the static structure factor, a collective 
observable amenable to direct observation in time-of-flight ex¬ 
periments for quantum spin chains realized by atom-optical 
quantum simulators. Finally, we have also briefly discussed 
the possibility of the direct observation of frustration, as quan¬ 
tified by our global-to-local incompatibility, by measuring the 
visibility in experiments with atomic interferometers. 

The investigation was carried out for one-dimensional sys¬ 
tems with competing nearest neighbor (NN) and next-to- 
nearest neighbor (NNN) interactions, but the methodological 
framework is general and in principle can be applied also to 
the study of systems with long-range interactions and/or to 
higher dimensions and different lattice geometries. 

The concepts introduced in the present work might be par¬ 
ticularly relevant to the investigation of systems that feature 
maximally quantum (nonlocal) and global orders, such as 
quantum spin liquids lt55il featuring topological order lf56tl . 
Topologically ordered phases do not break any symmetry of 
the many-body Hamiltonian, do not admit local order param¬ 
eters, and feature a robust ground-state topological degener¬ 
acy and a complex pattern of long-range ground-state entan¬ 
glement. A typical global signature of topological order is 


a nonvanishing sub-leading constant contribution to the area- 
law leading scaling of the ground-state block entanglement, as 
measured by the von Neumann entropy of the block reduced 
density matrix, known as topological entropy. 

On intuitive grounds, the universal measure of frustration, 
Eq. ©, should be perfectly suited for the characterization of 
topological order and other maximally quantum orders, as it 
measures precisely the degree of their incompatibility with the 
local order induced by the local interaction terms. 

The global-to-local incompatibility overlap is defined in 
any spatial dimension (at variance with topological entropy, 
which is not defined for one-dimensional systems), it is in 
principle computable in a clearly and unambiguously defined 
way (topological entropy has to be computed by complex and 
extremely delicate subtraction schemes), and is amenable to 
direct observation by strategies relying on time-of-flight or 
interferometric experiments (Von Neumann topological en¬ 
tropy is not directly accessible to experimental observation, 
and one should rely instead on the measurement of the two- 
Renyi topological entropy). On this basis, we are currently 
investigating the scaling behavior of frustration as a charac¬ 
terization of topological order JH. 

In general terms, quantum ground states can be classified 
according to their distinct patterns of long-range entangle¬ 
ment. This provides the most basic categorization of quan¬ 
tum phases of matter, more fundamental than Landau’s sym¬ 
metry breaking paradigm. It would then be worth pursuing 
the investigation of the relation between frustration, topolog¬ 
ical order, and patterns of genuine multipartite entanglement, 
along the lines suggested by some recent studies 1157145911 . In¬ 
deed, the observed analogy between two-site dimerization and 
single-site factorization suggests that fc-site ground-state fac¬ 
torization might be a useful tool to understand hierarchies of 
quantum many-body orders in terms of complex patterns of 
multipartite ground-state entanglement. 

Acknowledgments: - SMG and FI are grateful to Marcello 
Dalmonte for valuable discussions. SMG and FI acknowledge 
the FP7 Cooperation STREP Project EQuaM - Emulators of 
Quantum Frustrated Magnetism (GA n. 323714). SMG and 
BCH acknowledge the Austrian Science Fund (FWF-P23627- 
N16). 


[1] G. Toulouse, Commun. Phys. 2, 115 (1977); J. Vannimenus and 
G. Toulouse, J. Phys. C 10 , L537 (1977). 

[2] J. Villain, J. Phys. C 10 , 1717 (1977); 

[3] S. Kirkpatrick, Phys. Rev. B 16 , 4630 (1977). 

[4] K. Binder and A. P. Young, Rev. Mod. Phys. 58, 801 (1986), 
and references therein. 

[5] M. Mezard, G. Parisi, and M. A. Virasoro (eds.), Spin Glass 
Theory and Beyond (World Scientific, Singapore, 1987), and 
references therein. 

[6] C. Lacroix, P. Mendels, and F. Mila (eds.), Introduction to Frus¬ 
trated Magnetism (Springer-Verlag, Heidelberg-Berlin, 2011), 
and references therein. 

[7] H. T. Diep (ed.), Frustrated spin systems (World Scientific, Sin¬ 
gapore, 2013), and references therein. 

[8] C. M. Dawson and M. A. Nielsen, Phys. Rev. A 69, 052316 


(2004). 

[9] A. Ferraro, A. Garcla-Saez, and A. Acin, Phys. Rev. A 76, 
052321 (2007). 

[10] S. M. Giampaolo, G. Adesso, and F. Illuminati, Phys. Rev. Lett. 
104 , 207202 (2010). 

[11] M. M. Wolf, F. Verstraete, and J. I. Cirac, Int. J. Quant. Inf. 01, 
465 (2003). 

[12] N. de Beaudrap, M. Ohliger, T. J. Osborne, and J. Eisert, Phys. 
Rev. Lett. 105 , 060504 (2010). 

[13] J. Schnack, Dalton Trans. 39 , 4677 (2010). 

[14] S. M. Giampaolo, G. Gualdi, A. Monras, and F. Illuminati, 
Phys. Rev. Lett. 107 , 260602 (201 1 ). 

[15] U. Marzolino, S. M. Giampaolo, and F. Illuminati, Phys. Rev. 
A 88, 020301(R) (2013). 

[16] S. Hill and W. K. Wootters, Phys. Rev. Lett. 78 , 5022 (1997). 



11 


[17] W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 

[18] V. Coffman, J. Kundu, and W. K. Wootters, Phys. Rev. A 61, 
052306 (2000). 

[19] T. J. Osborne and F. Verstraete, Phys. Rev. Lett. 96, 220503 
(2006). 

[20] K. Kim, M.-S. Chang, S. Korenblit, R. Islam, E. E. Edwards, 
J. K. Freericks, G.-D. Lin, L.-M. Duan, and C. Monroe, Nature 
465, 590 (2010). 

[21] A. Bermudez, J. Almeida, K. Ott, H. Kaufmann, S. Ulm, U. 
Poschinger, F. Schmidt-Kaler, A. Retzker, and M. B. Plenio, 
New J. Phys. 14, 093042 (2012). 

[22] S. Zippilli, M. Johanning, S. M. Giampaolo, Ch. Wunderlich, 
and F. Illuminati, Phys. Rev. A 89, 042308 (2014). 

[23] T. GraB and M. Lewenstein, EPJ Quantum Technology 1, 8 
(2014). 

[24] C. Senko, P. Richerme, J. Smith, A. Lee, I. Cohen, A. Retzker, 
and C. Monroe, Phys. Rev. X 5, 021026 (2015). 

[25] L.-M. Duan, E. Dernier, and M. D. Lukin, Phys. Rev. Lett. 91, 
090402 (2003). 

[26] J. Simon, W. S. Bakr, R. Ma, M. E. Tai, P. M. Preiss, and M. 
Greiner, Nature 472, 307 (2011). 

[27] J. Struck, M. Weinberg, C. Olschliiger, P. Windpassinger, J. Si- 
monet, K. Sengstock, R. Hoppner, P. Hauke, A. Eckardt, M. 
Lewenstein, and L. Mathey, Nature Physics 9, 738 (2013). 

[28] M. Lewenstein, A. Sanpera, and V. Ahufinger, Ultracold 
atoms in optical lattices: simulating quantum many-body sys¬ 
tems (Oxford University Press, Oxford, 2012), and references 
therein. 

[29] D. G. Angelakis, M. F. Santos, and S. Bose, Phys. Rev. A 76, 
031805(R) (2007). 

[30] M. J. Hartmann. F. G. S. L. Brandao, and M. B. Plenio, Phys. 
Rev. Lett. 99, 160501 (2007). 

[31] A. A. Houck, H. E. Tiireci, and J. Koch, Nature Physics 8, 292 

( 2012 ). 

[32] S. Zippilli, M. Grajcar, E. Il’ichev, and F. Illuminati, Phys. Rev. 
A 91, 022315 (2015). 

[33] G. Chaboussant, Y. Fagot-Revurat, M.-H. Julien, M. E. Hanson, 
C. Berthier, M. Horvatic, L. P. Levy, and O. Piovesana, Phys. 
Rev. Lett. 80, 2713 (1998). 

[34] K. R. Brown, R. J. Clark, and I. L. Chuang, Phys. Rev. Lett. 97, 
050504 (2006). 

[35] F. Gerbier, A. Widera, S. Foiling, O. Mandel, T. Gericke, and I. 
Bloch, Phys. Rev. Lett. 95, 050404 (2005). 

[36] A. Hoffmann and A. Pelster, Phys. Rev. A 79, 053623 (2009). 

[37] P. M. Preiss, R. Ma, M. E. Tai, A. Lukin, M. Rispoli, P. Zu- 
pancic, Y. Lahini, R. Islam, and M. Greiner, Science 347, 1229 
(2015). 


[38] K. Uyanik and S. Turgut, Phys. Rev. A 81, 032306 (2010). 

[39] M. Blasone, F. DelT Anno, S. De Siena, and F. Illuminati, Phys. 
Rev. A 77, 062304 (2008). 

[40] M. Koashi and A. Winter, Phys. Rev. A 69, 022309 (2004). 

[41] S. M. Giampaolo, G. Adesso, and F. Illuminati, Phys. Rev. Lett. 
100, 197201 (2008). 

[42] S. M. Giampaolo, G. Adesso, and F. Illuminati, Phys. Rev. B 
79, 224434 (2009). 

[43] S. M. Giampaolo, S. Montangero, F. DelT Anno, S. De Siena, 
and F. Illuminati, Phys. Rev. B 88, 125142 (2013). 

[44] For a review of the J\ — J 2 model see U. Schollwock, J. Richter, 
D. J. J. Farnell, and R. F. Bishop (eds.): Quantum Magnetism. 
Lecture Notes in Physics, vol. 645 (Springer, Berlin, 2004). 

[45] C. K. Majumdar and D. Ghosh, J. Math. Phys. 10, 1388 (1969); 
C. K. Majumdar, J. Phys. C: Solid State Phys. 3, 911 (1969). 

[46] P. W. Anderson, Mater. Res. Bull. 8, 153 (1973); P. W. Ander¬ 
son, Science 235, 1196 (1987). 

[47] S. A. Kivelson, D. S. Rokhsar, and J. P. Sethna, Phys. Rev. B 
35, 8865 (1987); D. S. Rokhsar and S. A. Kivelson, Phys. Rev. 
Lett. 61, 2376 (1988). 

[48] I. Affleck, T. Kennedy, E. H. Lieb, and H. Tasaki, Phys. Rev. 
Lett. 59, 799 (1987); I. Affleck, T. Kennedy, E. H. Lieb, and H. 
Tasaki, Comm. Math. Phys. 115, 477 (1988). 

[49] M. Blasone, F. DelT Anno, S. De Siena, S. M. Giampaolo, and 
F. Illuminati, J. Phys.: Conf. Ser. 174, 012064 (2009). 

[50] M. Blasone, F. DelT Anno, S. De Siena, S. M. Giampaolo, and 
F. Illuminati, Phys. Scr. T140, 014016 (2010). 

[51] M. Dalmonte, S. M. Giampaolo, G. Zonzo, and F. Illuminati, 
Detecting topological order by quantum frustration, to appear. 

[52] A. K. Ekert, C. M. Alves, D. K. L. Oi, M. Horodecki, P. 
Horodecki, and L. C. Kwek, Phys. Rev. Lett. 88, 217901 
( 2002 ). 

[53] D. A. Abanin and E. Dernier, Phys. Rev. Lett. 109, 020504 
( 2012 ). 

[54] A. J. Daley, H. Pichler, J. Schachenmayer, and P. Zoller, Phys. 
Rev. Lett. 109, 020505 (2012). 

[55] L. Balents, Nature 464, 199 (2010), and references therein. 

[56] L. D. Carr (ed.). Understanding Quantum Phase Transitions 
(CRC Press, Taylor & Francis Group, Boca Raton, USA, 2011), 
and references therein. 

[57] S. M. Giampaolo and B. C. Hiesmayr, Phys. Rev. A 88, 052305 
(2013). 

[58] S. M. Giampaolo and B. C. Hiesmayr, New J. Phys. 16, 093033 
(2014). 

[59] L. Jindal, A. D. Rane, H. S. Dhar, A. Sen De, and U. Sen, Phys. 
Rev. A 89, 012316(2014). 


